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ABSTRACT 

I , Context. Be stars are fast-rotating early-type emission line stars. It is generally assumed that observed emission is generated in a 

. rotating disk-like envelope, as supported by the observed correlation between the stellar projected rotational velocity vsin; and the 

' width of the emission lines. Then, high-resolution spectroscopic observations of Balmer lines profiles play an important role in putting 

' constraints on Be stars modeling. 

Aims. We present Balmer lines spectroscopy for a sample of 48 Be stars. For most of them. Ha and Hfi have been observed more than 
two times, in a total period spanning almost two years between 2008 and 2009. 

Methods. Spectral synthesis of the Ha profile was performed following two steps: photospheric contribution was computed by using 
Kurucz's code ATLAS9 and SYNTHE, and disk emission was derived by the approach of Hummel & Vrancken (2000, A&A, 302, 
751). 

Results. For 26 out of 48 stars, a modeling of the total Ho- emission, i.e. photospheric absorption plus disk net emission, has been 
attempted. By this modeling we derived an estimation of the disk radius, as well as the inclination angle between the rotational axis 
with line of sight and the base density at the stellar equator. For the stars observed more than once, we also discuss the variability of 
Ha and H/3 for what concerns both the equivalent width and the spectral profile. We found 16 stars with variable equivalent width and 
7 stars with clear signs of profile variations. 
(~| ' Conclusions. For all the stars in our sample, we derive all the fundamental astrophysical quantities, such as, effective temperature, 

gravity, and projected rotational velocity. We found 13 stars whose equivalent width is variable with a confidence level greater than 
I I 80% and 7 object for which spectral profiles show change with time. According to the commonly used classification scheme, we 

O ■ classified 16 stars as belonging to class 1, 13 to class 2, 11 are shell stars, 6 objects do not show net emission, and 2 stars display 

transitions from class 1 and 2. For the class 1 stars, we confirm the correlation between v sin ; and peak separation. Concerning the 
C/3 ' geometry of the disk, we derived the base density at the stellar equator, the radius, and the inclination angle between rotational axis 

and line of sight. The maximum concentration of stars occurs for disk dimensions ranging in the interval of 6 to 8 stellar radii and for 
inclination angles going from 23° to 35°. 

Key words. Stars: emission-line. Be - stars: fundamental parameters 

oo 

\^ 1. Introduction estimated that Be stars rotate at a equatorial velocity equal to 80 

\^ % of the critical rotation velocity: 

• Classical Be stars are early-type B-type stars whose spectra have 

' one or more emission lines in the Balmer series. In particular, the =08 / ^'^* 

Ha emission line is typically the dominant feature in the spectra V R* 

_j of such stars, and many authors have modeled Ha line profiles t-, , ■ .... ■ t u uu 

, , , , , , The observed emission lines take a variety or shapes, which 

• • to understand the Be star phenomenon better. e u ■ u j u u u-i /Tmfei 

^ ■ ^ following the scheme proposed by Hanuschik (119961 1. range 

The emission lines observed in the spectra are explained in from wine bottle profiles, singly or double-peaked profiles, to 

terms of the recombination that occurs in a flattened circumstel- shell spectra, when the centi-al absorption must extend below the 

5^ ■ lar disk, according to the widely accepted model first proposed stellar- continuum flux. The various shapes are explained as a 

by Struve (|1931i l. The disk is a decretion disk; i.e., the source of dependency of /, the inclination angle of the star's rotation axis 

the disk material is the central star, generated by the equatorial to the observer's line of sight. In particular, shell profiles occur 

flow of stellar material. One of the key factors in creating the only when the disk is viewed equator-on (/ = 90°), while the 

disk is supposed to be the very high value of rotational velocity, single peak and wine botde occur only for near pole-on (/ = 0°) 

In fact. Be stars are known to have higher rotational velocities viewings, and double-peaked profiles occur at mid-inclination 

than a sample of normal B-type stars. From statistical consider- angles. 

ations on the v sin / distribution among Be stars, Porter (,1996j Double-peaked profiles have been observed both symmetric, 

which are the two peaks have the same intensity, and asymmet- 
ric, the peaks have different height over the continuum level. The 
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* Table B.l and all observed spectra are available in electronic form cuiTent theory is that asymmetry arises from one-aimed den- 

at the CDS via anonymous ftp to cdsai-c.u-strasbg.fr (130.79.128.5) or sity waves in the circumstellar disk, which is also known as 

via http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/ the global disk oscillation model. In this model, a one-armed 

** I wish to dedicate this paper to my child never bom, to keep track oscillation mode is superposed on an unperturbed, axisymmet- 

of his passage through my life. ric disk (Okazaki, 1997i) . Another aspect of Be stars emission 
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Table 1. Spectral types, luminosity classes, Stromgren pho- 
tometry, and derived effective temperatures for our program 
stars. Spectral types and luminosity classes are taken from the 
SIMBAD database 
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Fig. 1. Histogram showing the distribution of program stars as a 
function of spectral type. 



is their variability. For example, about one third of all double- 
peaked profiles exhibit changing asymmetry, with the so-called 
violet-to-red ratio (V/R) being cyclically variable on timescales 
of years to decades. 

Observations in different spectral regions help astronomers 
probe different regions of the stellar disk and then put constraints 
on modeling these stars. For example, recently Meilland et al. 
( l2007i ) have used VLTI/AMBER to observe a Arae in the Br^ 
line, which was constrained very strongly the rotational property 
of its disk, concluded that its rotation is purely Keplerian. 

In this paper we present an atlas of observed Ha and HyS 
spectral lines in a sample of bright Be stars. We payed particular 
attention to modeling the Ha profile and to the time variability 
of B aimer lines profiles. 

2. Observation and data reduction 

The idea that underlies this catalog is to create a homogeneous 
data set of stars observed with the same instrument, the telescope 
of the M.G. Fmcastoro station of INAF-Catania Astrophysical 
Observatory. For this purpose, we queried the "Catalogue of Be 
stars" compiled by Jaschek & Egret (1982 i selecting all the ob- 
jects with V<7 and observable at the latitude of the observatory, 
which means all the stars with 5 > -22° . The result of this query 
is the sample of 48 Be stars reported in Table [T] The limiting 
magnitude was chosen to obtain a good compromise between 
the exposure time and the signal-to-noise ratio. 

The present catalog is based on new spectroscopic obser- 
vations of all the stars in our sample, which spectral type are 
distributed between B 1 and AO (according to the histograms dis- 
played in Fig.[T]), and luminosities classes are III, IV, and V, as in 
the SIMBAD databas43- Some of these stars' spectra have never 
been published in other catalogs similar to ours, at least to our 
knowledge.late 

All the spectra of our program stars have been acquired 
with the 91 -cm telescope and FRESCO, the fiber-fed REOSC 
echelle spectrograph that allows spectra to be obtained in the 
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range of 4300-6800 A with a resolution R= 21 000. The spectra 
were recorded on a thinned, back-illuminated (SITE) CCD with 
1024x1024 pixels of 24 yum size, whose typical readout noise is 
of about 8 e" and the gain is 2.5 e"/ADU. All the spectra have 
been acquired during several observing runs spanning two years 
between 2008 and 2009. 

The reduction of spectra, which included the subtraction of 
the bias frame, trimming, correcting for the flat-field and the 
scattered light, the extraction of the orders, and the wavelength 
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Table 2. Adopted astrophysical quantities for our stars and final best-fit parameters. For each star we report, the number of spectra 
collected, their classification (abs means absorption lines only), the projected equatorial rotational velocity as measured by us in 
our spectra, mass and radius as estimated from Drilling & Landolt (119991 ), peaks separation, critical velocity, equatorial velocity, 
inclination angle, base density at the stellar equator, and the estimated disk radius. In the last column we also reported the references 
of other similar catalogs in which spectra of the star have been published: 1) Hanuschik ( 119881 ), 2) Hanuschik et al. ( 1 19961 1, 3) 
Hanuschik ( [19961 ), 4) Slettebak et al. ( [1992] ), 5) Saad et al. ([20061 ), 6) Silaj et al. ([20T0t . 
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calibration, was done by using the NOAO/IRAF packagfl The 
amount of scattered light correction was about 10 ADU. After 
dividing the extracted spectra by flat-field, the residual shape 
the spectrum was removed by dividing each spectral order by 
a Legendre function of a low order. Typical S/N of our spectra is 



^ IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Universities 
for Research in Astronomy, Inc. 



~ 100. For some stars this limit has not been reached because its 
apparent magnitude is close to Va!7, in which case the S/N was 
about 50. 

Finally, the IRAF package rvcorrect was used to include the 
velocity correction due to the Earth's motion, which moved the 
spectra into the heliocentric rest frame. The task splot and its 
facilities were used to measure the peaks separation in the Ha 
profile. Errors in the pixels position were converted in errors on 
the separations, and were evaluated in w 20 km s"' . 
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Table 3. Program stars for which we have more than two spec- 
tra. We reported for both lines mean equivalent widths, their F 
values, computed according Eq. |3] and the confidence levels of 
the detected variabihty. 
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For each spectrum we measured the equivalent widths (in- 
cluding underlying absorption) of both Ha and Hy6, where a neg- 
ative value means the corresponding line shows net emission. All 
the measured equivalent widths are reported in Table IB. fl 



3. Classification and fit of the t-io' line profiles 

Almost all the stars in our sample were found with emission in 
the Ha. Then, considering the shape of their profile and accord- 
ing to the classification scheme proposed by Hanuschik (I1988I I. 
we classified our stars as belonging to 

- class 1 when they exhibit a rather symmetrical double peak 
structure with V/R x 1, 

- class 2 when they have an asymmetric single peak or a dom- 
inant peak with a much weaker secondary peaked, 

- shell stars when the central reverse is deeper than the contin- 
uum level, 

- abs when there is not emission above the continuum level. 

Our classification is reported in the fourth column of Table |2] 

Only for the stars belonging to class 1, plus stars from other 
classes but with V/R ratio close to unity, we attempt an estima- 
tion of the disk dimension. The approach we used was to min- 
imize the difference between observed and synthetic profiles. 



computed in two separate steps. First of all, we calculated the 
photospheric Ha, then the contribution due to the net emission of 
the disk and then we added these two synthetic profiles obtained 
separately. These two steps are described in the following: 

Computation of the photospheric profile 
We first computed the photospheric Ha profiles for all our 
program stars. They were generated in three steps: i) first, 
we computed an LTE model atmosphere using the ATLAS9 
code (Kurucz 1993), ii) the stellar spectrum was then syn- 
thesized using SYNTHE (Kurucz & Avrett [19811) . and Hi) 
the spectrum was convolved with the instrumental and rota- 
tional profiles. 

First of all we had to obtain an estimation of the effective 
temperature for each target. Considering that the continuum 
energy distribution of Be stars is typical of normal early-type 
stars both in the visual and UV, but not in the IR, where 
an excess could be present because of the hot circumstellar 
dust (Zickgraf 2000), effective temperatures were computed 
from Stromgren photometry (Hauck & Mermilliodil998 1 us- 
ing the algorithm coded by Moon (1 19851 1. with the exception 
of eight stars for which photometry is not available. This 
method is allowed because it does not involve any IR fil- 
ter For seven of them we adopted the temperatures from 
the literature: HD 10516, HD 37202, and HD 44458 from 
Soubiran et al. (2010), HD 58020 and HD 71072 from Hohle 
et al. (120101) . HD 162428 from Moujtahid et al. (1999), and 
HD 212571 from Wu et al. (2011) . while for HD 53416 
we derived an estimation of temperature from spectral 
type and the calibration by Kenyon & Hartman ( |19951 l. 
Since our targets have luminosity class IV/V (as reported 
in the SIMBAD database), we fixed the surface gravity 
to log g = 4.0, except for HD 11415, HD 37490, HD 45542, 
HD 50658, HD 109387, HD 19391 1, and HD 217050 (lumi- 
nosity class III) for which log ^ = 3.0 has been preferred. 
Radii and masses were adopted following the calibration in 
Drilling & Landolt (119991 ). Assuming these atmospheric pa- 
rameters, we computed the vsin / of each star by spectral 
synthesis of the observed Mgi /14481 A. This line was cho- 
sen because in the spectral range of our targets, it reaches its 
maximum depth and therefore it is better suited to determin- 
ing the rotational velocity. Errors on the projected rotational 
velocities are w 15 km s"'. 
Computation of net disk emission 

We have adopted the Be disk model approach of Hummel 
& Vrancken ( 2000) that is based on models developed by 
Home & Marsh (119861 ) and Home (119951 ) for accretion disks 
in cataclysmic variables. The disk is assumed to be axisym- 
metric and centered over the equator of the underlying star, 
and the gas density varies as 



p(R,Z) = pQR'"exp 



2 \mR)l 



where R and Z are the radial and vertical cylindrical coordi- 
nates (in units of stellar radii), po is the base density at the 
stellar equator, n a radial density exponent, and H(R) the disk 
vertical scale height. The neutral hydrogen population within 
the disk is found by equating the photo-ionization and re- 
combination rates (Gies et al. 2007 ). The disk gas is assumed 
to be isothermal and related to the stellar effective tempera- 
ture Teff by T^ = 0.6 Teff (Carciofi & Bjorkman 2006 ). 
This approach take the contribution of the central star's finite 
size on the Ha line formation process into account, i.e. the 
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obscuration of the disk by the central star at any given incli- 
nation. The numerical model represents the disk by a large 
grid of azimuthal and radial surface elements, and the equa- 
tion of transfer is solved along a ray through the center of 
each element according to 



■') + e-^' 



where I,i is the derived specific intensity, the source func- 
tion for the disk gas (taken as the Planck function for the 
disk temperature Tj), the specific intensity for the Ha of 
the star, and the integrated optical depth along the ray. The 
first term applies to all the disk area elements that are unoc- 
culted by the star, while the second term applies to all ele- 
ments that correspond to the projected photospheric disk of 
the star. The absorption line adopted in is Doppler-shifted 
according to solid-body rotation for the photospheric posi- 
tion in a star that is rotating at 80% of the critical value. 
Electron scattering is not taken into account in the line pro- 
file computation. So we do not expect to reproduce the wings 
of strong lines well, since for these lines the broadening of 
the wings due to the electron scattering can not be neglected. 
Disk kinematic is taken into account using a rotational ve- 
locity law written as 

v(R) = v:^,R-j 

(Hutchings 1 1 9701 ). where R represents the radial coordinate 
that has its origin at the center of the star, and Vj^, denotes 
the actual rotational velocity at the stellar surface. The expo- 
nent ranges from j = 1/2 for pure Keplerian rotation and 
7 = 1 cotTesponding to conservation of angular momentum. 
Likewise the value of j is still matter of debate, recent stud- 
ies seem to converge toward the Keplerian value (Hummel 
& Vrancken 120001 Meilland et al. 2007 ). Thus, in this study 
we assumed the disk to be in pure Keplerian rotation. 

Once we obtained and combined these two contributions, we 
started the minimization algorithm using as goodness-of-fit test 
the parameter 



X 



where is the total number of points, /obs and /th are the inten- 
sities of the observed and computed profiles, respectively, and 
5/obs is the photon noise. 

As initial guesses for the inclination angle ; and for the disk 
radius, we used the equations 



Veq sin i = 0.8 



sm I 



where Veq sin / is the value measured in our spectra, and 

AVpeak 

2 V sm i " 



(1) 



(2) 



(equation [6] in Hanuschik et al. 09881 1) where Avpeak is the sep- 
aration between violet and red peaks, as measured in each profile 
and reported in Table |2] 

Then, fine tuning was carried out using the amoeba mini- 
mizatior0 algorithm between observed and computed profiles. 



^ The amoeba routine implements the simplex method of Nelder & 
Mead jl965D . 




(km s-> 



Fig. 2. CotTelation between measured v sin / and peaks separa- 
tion. The points have been grouped on the basis of the mean 
equivalent width: EW<-20 A (filled squares), -20A<EW<- 
15 A (filled circles), -10A<EW<-5A (filled triangles), - 
5 A < EW < A (open circles), and EW > A (open squares). 



In our procedure two assumptions have been made: radial den- 
sity exponent has been fixed to n-3 and, as stated before, the 
Keplerian rotation of the disk has been considered (j=l/2). 
The first hypothesis, regarding the value of the density expo- 
nent can be justified by considering the work of Grundstrom & 
Gies (2006). These authors computed several theoretical curves 
that described the dimension of the disk radius as a function of 
the Ha equivalent width, for different values of the inclination 
angle i and different values of n. They concluded that the over- 
all shape of those curves for different n and equal i are almost 
the same, since it is a small difference of a; 3% in correspon- 
dence of equivalent widths between -2 and -15 A when n change 
from 3 to 3.5. They then suggest that the particular choice of n 
is not as important as the choice of the right i. Moreover, Porter 
& Rivinius ( I2003I I from IR flux excess in Be stars suggested that 
n falls in the interval n - 2-^4. Thus, on the basis of these 
results, we fixed the value of the density exponent to the mid- 
dle value of n = 3. Total Ha profiles, star-ndisk, are presented in 
Figs.|Al]|A21 and|A3] 

To derive an estimation of the disk radius, we used the 
method developed by Grundstrom & Gies (2006 1 to form a syn- 
thetic image of the system star+disk in the plane of sky by sum- 
ming the intensity over a 2.8 nm band centered on Ha. We 
collapsed this image along the projected major axis to get the 
summed spatial intensity, and we adopted the value for which 
the summed intensity drops to half its maximum value as effec- 
tive disk radius. 

We estimated the erTors on the disk dimensions and on the 
inclination angles to be +2 R, and +3°, respectively. These de- 
terminations have been estimated by varying in Eqs. [T]and in|2] 
the observed quantities v sin / and peak separation by their ex- 
perimental errors and considering as uncertainties the semi- 
amplitude of this variation. 

All the adopted and derived parameters are reported in 
Table|2] 
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4. Hof and H/3 Variability 

Usually Be stars display variability in their equivalent width 
(EW) and/or in their spectral profile. 

To find whether a star presents equivalent width variation, 
we applied to both Ha and H/3 equivalent widths the statisti- 
cal method called F-test. When more than one observation was 
present for a given target, we calculated the amplitude of the 
variation, AEW, and the standard deviation of the sample using 

where N is the number of points, and the (N-1) corresponds to 
the degree of freedom used for the F-test. Having obtained AEW 
and cr, we calculated a simple observable to assess the variability 
for a given target using a ratio of the form: 

F . (3) 
cr 

This simple ratio represents the number of times that the ampli- 
tude of the variation is greater than the standard deviation. 

To determine whether or not this number is meaningful and 
whether the star shows variability, we evaluated the correspond- 
ing confidence level; see last column of Table [3] We considered 
all the stars for which C > 80% as definitively variable (13 stars, 
that is 41% of the sample), while we cannot say anything for the 
target with C < 50% (22% of the sample) due to the small num- 
ber of points. The 37% of the stars, are probably variable but 
within a confidence level ranging from 50% to 80%. In any case, 
all the targets show the same behavior in both spectral lines. 

In this atlas we show the full set of our profiles (Ha and 
Hj3) for the 48 observed stars. In Figs. IB.ll to lB3] we show the 
Balmer profiles observed for program stars with multiple mea- 
surements. In each panel we reported for each profile also the 
last four digits of the heliocentric Julian date of the observation. 
In some cases to improve the visualization, profiles have been 
blown up according to 

F',^[(F,-l)*k] + l (4) 

where denotes the observed flux, the displayed flux, and k 
is the magnification factor. 

In Fig. lB.7l we show all the stars for which we collected one 
spectrum only. Most of the stars have been measured several 
times during the 2008/2009 period, but only few stars showed 
significant variations in their spectral profiles. These objects are 
discussed separately in the following: 

HD 37202 - This star has been observed in two nights sepa- 
rated by 205 days. It is clearly visible an increase in flux in the 
violet peak in both spectral lines, even if it is more evident in 

HD 41335 - As for the previous object, only two observations 
have been acquired for this star in a range of 207 days. The star 
exhibits V/R variations in both lines. 

HD 58050 - The tiiple-peaked structure visible in the Ha ob- 
served in the first night is missing in the other two spectra. No 
features are seen in the H/3. 

HD 109387 - This star was observed for 13 nights in a pe- 
riod spanning 110 days. The top part of the Ha emission profiles 
shows irregular variability, with back and forth changing from 
class 2 toward class 1, but no evident sign of day -by-day vari- 
ability has been observed. No variations have been detected in 
Hy6 double-peaked profile. 



HD 142926 - During 110 days this shell star was observed 
10 times showing a slight V/R variability in the Ha. 

HD 143275 - In the four spectra acquired by us, S Sco shows 
important changes in the emission of Ha. The V/R changes over 
the observational period and it shows a flat core in the first spec- 
trum. 

HD 164284 - Observed for six times in 536 days, this star ex- 
hibits equivalent width variability. The double-peaked emission 
of the Ha profile decrease with time, although the V/R remains 
constant and w 1 . Also the H/3 shows a change in the shape, since 
it is the last profile without any peak. 

HD 183362 - This star shows an increased emission level on 
the red side of its profile. 

HD 183656 - V/R variability has been observed, both in Ha 
and in Hp profiles, in the six spectra acquired in a range of 358 
days. 

HD 187567 - This star show variability in the Ha profile, and 
evolves from class 1 toward class 2. 

HD 189687 - In the Ha, this star does not show any sign of 
variability for the first month of observations. In the last spec- 
trum taken after 32 days from the second to last, it starts to show 
an increase of the flux in the red peak. No H/3 variability has 
been detected. 

HD 191610 - This star shows a gradual increment of the flux 
in the red peak. 

5. Discussion and conclusions 

In this paper we presented a homogeneous sample of Ha and 
H/3 line profiles observed in 32 Be stars, which show emission at 
least in the Ha line. According to Hanuschik ( 1988), we classi- 
fied our targets on the basis of the following scheme: 

- 16 stars, 33% of the sample, in class 1; 

- 13 stars, 27% of the sample, in class 2; 

- 1 1 stars, 23% of the total, has been classified as shell stars; 

- 6 stars, 12% of the total, do not show net emission. 

In this list we do not include the two stars that show any phase 
transition between classes 1 and 2. 

This frequency distribution shows that the majority of our 
sample of 48 Be stars, randomly distributed in spectral type, be- 
longs to class 1 profiles. Regarding the 13 stars classified as class 
2, seven of them are single peak, while six show structured pro- 
files. 

Two stars showed variability from one class to another. 
HD 187567 has undergone an evolution from class 1 to class 2, 
while the behavior of HD 109387 is more complicated. In 110 
days, this star has been observed 13 times, and it showed a tran- 
sition from class 2 to class 1 and back again to class 2. 

In Fig. [2] we compared the behavior of the Ha peaks sep- 
aration to V sin / for the stars of the class 1 (included the class 
2 HD 60855). A linear correlation seems to exist, as expected 
from the work of Hanuschik et al. (il988) . although two stars 
discarded from this trend, namely HD 191610 and HD 212571. 
To verify this correlation we computed the Pearson r coefficient, 
obtaining r = 0.72 and the linear fit given by the equation: 

y = (0.56 + 0.13) ■ X + (29.94 + 26.80). (5) 

Thus, we confirm that vsin/ and Ha peaks separation are 
in linear correlation, confirming the disk-like geometry of Be 
star envelopes and, probably this assumption is not valid for 
HD 191610 and HD212571. 
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Fig. 3. Histograms of distribution of stars as a function of the dimension of their disks (left panel) and of the inclination angles (right 
panel). To build these histograms we chose a bin size of 2 R, and 3°, respectively. 



For all the stars belonging to the class 1, we attempted to 
model the emission with the purpose of deriving some param- 
eters such inclination angle, base density at the stellar equator 
and disk radius. In Fig.|3]we show the distribution of stars as a 
function of the disk radius (right panel) and of the inclination an- 
gle (left panel). The histograms were built considering a binning 
equal to the estimated eiTors, that is, 2R, and 3° in the disk di- 
mension and inclination angle, respectively. They show as there 
is a major concentration of stars for disk around 6 H- 8 R, (about 
17% of our sample) and for angles around 23° -H 35° (about 28% 
of the sample). 

Moreover, with the aim of inferring line profile variabil- 
ity, for most of the stars of our sample we obtained more than 
one spectrum in a period spanning two years between 2008 and 
2009. All but seven stars, those discussed in Sect.lH do not show 
any evident sign of variability in both Balmer lines. 
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G. Catanzaro: Spectroscopic atlas of Hor and US in a sample of northern Be stars 



Table B.l. Ha and equivalent width as measured in our spectra. Negative values mean that the lines show net emission. 
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Fig. B.l. Line profiles of Hor and HjS observed for our program stars. Spectra have been shifted along the vertical axis according to 
^e Julian Day of their observation, whose last four digits are reported. In parenthesis we report the magnification factors apphed to 
each profiles in order to improve the visualization. 
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Fig.B.2. As inFig.lO 
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Fig.B.3. As inFig.lBH 
16 



G. Catanzaro: Spectroscopic atlas of Ha and H/3 in a sample of northern Be stars 



HD164447 HD183362 

JD = 2450000 + JD = 2450000 + 




-10 -5 5 10 -10 -5 5 10 15 -10 -5 5 10 -10 -5 5 10 15 
Velocity (*100 km/s) Velocity (*100 km/s) Velocity (*100 km/s) Velocity (*100 km/s) 

HD183656 HD183914 

JD = 2450000 + JD = 2450000 + 




-10 -5 5 10 -10 -5 5 10 15 -10 -5 5 10 -10 -5 5 10 15 
Velocity (*100 km/s) Velocity (*100 km/s) Velocity (*100 km/s) Velocity (*100 km/s) 

HD187567 HD189687 

JD = 2450000 + JD = 2450000 + 




-10 -5 5 10 -10 -5 5 10 15 -10 -5 5 10 -10 -5 5 10 15 

Velocity (*100 km/s) Velocity (*100 km/s) Velocity (*100 km/s) Velocity (*100 km/s) 



Fig.B.4. As inFig.lO 
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Fig.B.6. As inFig.lBH 
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Fig. B.7. Ha and H/3 profile for stars for which we only have one spectrum. In the boxes we display shell and class 1 stars, class 2 
stars, and stars showing absorption, respectively. 
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